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Abstract 

We show that a unified approach to the perturbative evolution of structure functions which 
sums all logarithms of Q 2 and 1/x at leading and next-to-leading order yields results in 
full agreement with the 1993 HERA data for F2. This makes it possible to determine a s 
surprisingly accurately from these data alone. 
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At small x the structure function F^x.Q 2 ) measured recently at HERA de- 
pends in effect on two scales, Q 2 and s = ("4r)Q 2 , instead of just the single scale Q 2 
relevant at large x. The renormalization group equations which govern the perturba- 
tive evolution of singlet parton distribution functions thus describe evolution in both s 
and Q 2 (or equivalently 1/x ~ s/Q 2 and t = lnQ 2 /A 2 ), and should sum up all leading 
(and subleading) logarithms of both Q 2 and s (or Q 2 and 1/x) . In fact it is possible to 
construct unified perturbative evolution equations for leading twist singlet parton distribu- 
tion functions which are straightforward generalizations of the usual singlet Altarelli-Parisi 
equations, but with new splitting functions which explicitly incorporate all the logarithms 
of 1/x. The resulting equations^] may be used to calculate perturbative evolution down 
to arbitrarily small values of x, provided only that Q 2 is sufficiently large that higher twist 
effects may be ignored. All nonperturbative effects (at both large and small x) are then 
factorized into the initial condition at to. 

The singlet splitting functions in the the double leading expansion scheme, which 
treats the two scales symmetrically (see ref.Q), take the form P i: >(x,t) = P^ (x) + 
^^LoOM) + where if G (x) = if (x) + P*i(x, t), P^ LO = if (x) + P% (x, t), 
If (x) and If (x) are the usual one and two loop splitting functions, while If (x,t) and 
Pgi(x,t) are (convergent) sums of leading and subleading singularities respectively: for 
x > xq they are all zero, while for x < xq P§ 9 (x,i) = P]f q (x,t) = 0, but 

00 

xP^(x,t) = §LxP°*(x,t) =2C A J2an^ T yM t ) n ~ 1 C-\ 

n=4 

00 

zPg(z,t) = ^xp™(x,t) = §r a n / £a n ^A.(t) n - 1 r- 2 , (i) 

Tl = 3 

00 

xPg(x,t) = %L X Pg(x,t) =2C A J2bn J ^ W \s(tr- 1 C- 2 , 

n=3 

where X s (t) = 41n2^a s (t), a s (t) is the two-loop running coupling, and £ = log 2a . 

The coefficients (a n , a n ) may be found in ref.@: they are computed using expressions 
derived in ref.[|J. The coefficients b n (and the color-charge relation between Pf s 9 and Pf/) 
can be fixed uniquely at NLO by requiring momentum conservation: b n = —a n — -^P- a n . 
This fixes the arbitrariness implicit in alternative schemes in which momentum is conserved 
separately for quarks and gluons^, and furthermore makes the explicit computation of 
P§jj and Pf s 9 unnecessary at NLO: subleading corrections to the BFKL kernel would only 
be necessary for NNLO computations. The parameter xq, which marks the boundary 



between large and small x regions, must eventually be fixed empirically: we vary it in the 
range < xo < 0.1, considering larger values of xo to be unreasonable. 

Each of the series for the splitting functions is uniformly convergent on any closed 
interval of x which excludes the origin, and can thus be used to extend perturbative 
evolution down to arbitrarily small values of a?Q. It follows that their Mellin transforms 
(the anomalous dimensions 77V ) have no physically relevant singularities beyond those 
present at fixed order in a s (so in particular any cuts generated by all order resummations 
may in practice be ignored). 

The asymptotic behaviour of the parton distribution functions, as solutions to these 
unified evolution equations, depends on the particular limit adopted in the x-t plane. In the 
Bjorken limit t — > 00, with x fixed, the most important parts of the splitting functions are 
Pi to LO and P2 at NLO. The dependence on x is essentially nonperturbative, being given 
by the initial condition at to, while the dependence on t is perturbative. In the Regge limit 
£ — ► 00 at fixed t, the most important parts of the splitting functions are the leading and 
subleading singularities P s and P ss . In this limit F% eventually rises as x~ x ^\ the power 
X(t) being given by Regge theory for t ~ to, but calculable perturbatively for t sufficiently 
large (X(t) = A s (t))||. Finally, in the double scaling limit £ — > 00, £ = In ^ — > 00, 
rising as some power of £, but slower than an exponential, the asymptotic form of singlet 
parton distribution functions may be determined completely in perturbation theory || J7| [Q. 
In particular the structure function Fz{x, t) then grows as Rp 1 = A/'(j"2 / o~ 1 e 27<J - <5c , where 
a = y/£C, p = \/|7C) 7 = f an d 8 = fg- This leads to a double asymptotic scaling in a 
and p which is indeed seen experimentally 0@(see fig. 1). In the double scaling limit (and 
thus in the kinematic regime studied at HERA) the most important parts of the splitting 
functions are thus the pivotal leading singularity of Pi (which gives 7) and the subleading 
term (which gives 5) at LO, and the leading singularity of P2 at NLO. 

Two loop calculations of F^ in the double scaling region were performed in ref.||: 
here we give similar calculations in the double leading scheme. The solution of the NLO 
evolution equations is straightforward, although care must be taken to linearize all sub- 
leading corrections, to avoid spurious sub-subleading terms in the solution. As in [[| we 
use as input the MRS distribution DO' with a small-x tail parameterized by the single pa- 
rameter A(Qo) at the starting scale Qq. The two parameters X(Qo) and a s are then fitted 
by minimising the total x 2 of the evolved distribution to the HERA data. The resulting 
best- fit parameters are given in table 1, and the corresponding distributions displayed in 
figure 1, in the form of two scaling plots]?]]. The two loop curve (i.e. that with xq = 0) 
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Figure 1: Scaling plots of RfF% vs a and p. All 1993 HERA data passing 
the cuts a, p > 1, Q 2 < (2rrib) 2 are included, and Rf is thus evaluated with 
nj = 4. The diamonds and squares are ZEUS[|J and H1[Q] data, renormal- 
ized by 1.02 and 0.95 respectively. The curves are the simple double scaling 
prediction Q (dotted), and two NLO double leading calculations with ex- 
tremal values of xq\ xq = (i.e. two loops only) (solid) and xq = 0.1 
(dashed). The curves in the a-plot have p = 2.2: those on the p-plot 
a = 1.7. 



now fits the data significantly better than the scaling curve (which has x 2 = 71/121), 
essentially because the most singular part of the two loop correction reduces the slope of 
the a plot a little at lower values of a, a tendency which is evident in the data. The curve 
with xo = 0.1 is almost identical to that with xo = 0, differing only at large values of p; 
as explained in ref. 0, most of the effect of the higher order logarithmic singularities can 
be absorbed by suitably adjusting the boundary condition, either by changing Q at fixed 
A, or (as here) A at fixed Qq. If one were to keep the boundary condition fixed one might 
instead have the impression that the higher loop singularities have a large effect on in 
the HERA region. 

Whereas the higher order singularities only affect Fi close to the nonperturbative 
boundary at t = to, they can have a relatively large effect on the size of the quark and 
gluon distributions at small x and large t. To search for this effect (and, if it is found, to 
measure xq) it would be very interesting to have an independent measurement of the gluon 
distribution, either through Fl or heavy quark production. If xq is as large as 0.1, such a 
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102% 95% 


-0.23 ±0.03 


0.122 ±0.002 


58.4/120 


0.1 


99% 93% 


-0.06 ±0.06 


0.114 ±0.004 


64.9/120 



Table: Fitted parameters and % 2 for the two fitted curves in figure 1. Sta- 
tistical and systematic errors for each data have been added in quadrature, 
and the normalizations of the two experiments fitted within their stated 
uncertainties of ±3.5% and ±4.5% respectively (lj ||. 




measurement should find a significantly smaller gluon distribution than that expected at 
two loops. 

Since the rise in at small x and large Q 2 is being driven essentially by the triple 
gluon vertex, it depends rather strongly on the value of a s . It should thus be a good 
place to measure a s , requiring far lower statistics than at large x, where one has to search 
instead for small violations of Bjorken scaling|TtJ. Indeed the 1993 HERA data are already 
sufficient for such a determination, as is apparent from the two parameter fits of shown 
in the table (see also fig. 2), and the fact that the two loop correction can be seen in the 
data. A preliminary error analysis gives 



a s (M z ) = 0.120 ± 0.005(exp) ± 0.010(th) 



(2) 



a more detailed determination will be presented elsewhere. As at large x the theoretical 
error is dominated by the renormalization scale uncertainty, with an additional uncertainty 
here due to the unknown value of xq. Higher twist effects however seem to be very small. 
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